Shell model description of the ^^C dating [3 decay 
with Brown-Rho-scaled NN interactions 
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We present shell model calculations for the /3-decay of ^''C to the '^^N ground state, treating 
the states of the A = multiplet as two Op holes in an ^^O core. We employ low-momentum 
nucleon-nucleon (NN) interactions derived from the realistic Bonn-B potential and find that the 
Gamow- Teller (GT) matrix element is too large to describe the known lifetime. By using a modified 
version of this potential that incorporates the effects of Brown-Rho scaling medium modifications, 
we find that the GT matrix element vanishes for a nuclear density around 85% that of nuclear 
matter. We find that the splitting between the ( J", T) = (1+, 0) and (J", T) = (0+, 1) states in ^^N 
is improved using the medium-modified Bonn-B potential and that the transition strengths from 
excited states of ^'*C to the ^'*N ground state are compatible with recent experiments. 



The beta decay of ^^C to the ^''N ground state has long 
been recognized as a unique problem in nuclear structure. 
Its connection to the radiocarbon dating method, which 
has had a significant impact across many areas of science, 
makes the decay of broad interest even beyond nuclear 
physics. But a priori one would not expect the beta 
decay of ^"^C to be a good transition for radiocarbon dat- 
ing over archaeological times, because the quantum num- 
bers of the initial state [J^ ,T) = (0+, 1) and final state 
( J'^, T) = 0) satisfy the selection rules for an allowed 
Gamow- Teller transition. The expected half-life would 
therefore be on the order of hours, far from the unusually 
long value of 5730 years 1] observed in nature. The corre- 
sponding nuclear transition matrix element is very small 
(~ 2 X 10~^) and is expected to result from an accidental 
cancellation among the different components contribut- 
ing to the transition amplitude. This decay has therefore 
been used to investigate phenomena not normally con- 
sidered in studies of allowed transitions, such as meson 
exchange currents , relativistic effects and con- 
figuration mixing [5|,|6|. Of broader importance, however, 
is that this decay provides a very sensitive test for the 
in-medium nuclear interaction and in particular for the 
current efforts to extend the microscopic description of 
the nuclear force beyond that of a static two-body poten- 
tial fit to the experimental data on two-nucleon systems. 
One such approach is to include hadronic medium mod- 
ifications, in which the masses of mesons and nucleons 
are altered at finite density due to the partial restoration 
of chiral symmetry 0, H, Q or many-body interactions 
with either intermediate nucleon-antinucleon excitations 
[lot or resonance-hole excitations [lH. These effects are 
traditionally incorporated in models of the three-nucleon 
force (3NF), which have been well-tested in ah initio nu- 
clear structure calculations of light nuclei 12, 1^. 



In this Letter we suggest that a large part of the 
observed ^"^C beta decay suppression arises from in- 
medium modifications to the nuclear interaction. We 
study the problem from the perspective of Brown-Rho 



scaling (BRS) [li, [ist , which was the first model to make 
a comprehensive prediction for the masses of hadrons at 
finite density. In BRS the masses of nucleons and most 
light mesons (except the pion whose mass is protected 
by its Goldstone boson nature) decrease at finite density 
as the ratio of the in-medium to free-space pion decay 
constant: 
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where qa is the axial coupling constant, <I> is a func- 
tion of the nuclear density n with $(rio) — 0.8 at nu- 
clear matter density, and the star indicates in-medium 
values of the given quantities. Since all realistic models 
of the NN interaction are based on meson exchange and 
fit to only free-space data, eq. H]) prescribes how to con- 
struct a density-dependent nuclear interaction that ac- 
counts for hadronic medium modifications. This program 
has been carried out in several previous studies of sym- 
metric nuclear matter 1^ 17 1, where it was found that 



one could well describe saturation and several bulk equi- 
librium properties of nuclear matter using such Brown- 
Rho-scaled NN interactions. 

The case of the ^''C beta decay provides a nearly ideal 
situation in nuclear structure physics for testing the hy- 
pothesis of Brown-Rho scaling. Just below a double shell 
closure, the valence nucleons of ^^C inhabit a region with 
a large nuclear density. But more important is the sen- 
sitivity of this GT matrix element to the nuclear tensor 
force, which as articulated by Zamick and collaborators 
(isl ligj is one of the few instances in nuclear structure 
where the role of the tensor force is clearly revealed. In 
fact, with a residual interaction consisting of only cen- 
tral and spin-orbit forces it is not possible to achieve 
a vanishing matrix element in a pure configuration 
pot . Jancovici and Talmi 21 1 showed that by including 
a strong tensor force one could construct an interaction 
which reproduces the lifetime of ^''C as well as the mag- 
netic moment and electric quadrupole moment of ^''N, 
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although agreement with the known spectroscopic data 
was unsatisfactory. 

The most important contributions to the tensor force 
come from tt and p meson exchange, which act opposite 
to each other: 
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Since the p meson mass is expected to decrease substan- 
tially at nuclear matter density while the tt mass remains 
relatively constant, an unambiguous prediction of BRS is 
the decreasing of the tensor force at finite density, which 
should be clearly seen in the GT matrix element. In 
fact, recent shell model calculations 22] performed in a 
larger model space consisting of + 2huj excitations 
have shown that the /3-decay suppression requires the in- 
medium tensor force to be weaker and the in-medium 
spin-orbit force to be stronger in comparison to a typical 
G-matrix calculation starting with a realistic NN inter- 
action. We show in Fig. [T] the radial part of the tensor 
interaction V'^{r) = Vj{r) + V^{r) at zero density and 
nuclear matter density assuming that m* (jiq) /nip = 0.80. 
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FIG. 1: The radial part of the nuclear tensor force given in eq. 
([2]) from TT and p meson exchange at zero density and nuclear 
matter density under the assumption of BRS. 

Experiments to determine the properties of hadrons in 
medium have been performed for all of the light mesons 
important in nuclear structure physics. Studies of deeply- 
bound pionic atoms ^23,] find only a small increase in the 
7r~ mass at nuclear matter density and a related decrease 
in the n^ mass. Experimental information on the scalar 
and vector particles comes from mass distribution mea- 
surements of in-medium decay processes. Recent pho- 
toproduction experiments of correlated pions in the 
T = J — channel (cr meson) have found that the distri- 
bution is shifted to lower masses in medium. The vector 



mesons have been the most widely studied. Whereas the 
situation is clear with the w meson, the mass of which 
drops by ~ 14% at nuclear matter density [i^, with the 
p meson it is still unclear 2^ 23|. We believe that our 



present study tests the decrease in p mass more simply. 

Today there are a number of high precision NN interac- 
tions based solely on one-boson exchange. In the present 
work we use the Bonn-B potential [i^ which includes the 
exchange of the tt, rj, a, oq, p, and lu mesons. In [16| the 
consequences of BRS on the free-space NN interaction 
were incorporated into the Bonn-B potential and shown 
to reproduce the saturation properties of nuclear mat- 
ter in a Dirac-Brueckner-Hartree-Fock calculation. The 
masses of the pseudoscalar mesons were unchanged, and 
the vector meson masses as well as the corresponding 
form factor cutoffs were decreased according to 



A* n 
= — = 1 - 0.15 — . 
A no 



(3) 



The medium-modified (MM) Bonn-B potential is unique 
in its microscopic treatment of the scalar a parti- 
cle as correlated 27r exchange. Finite density effects 
arise through medium modifications to the exchanged p 
mesons in the pionic s-wave interaction as well as through 
the dressing of the in-medium pion propagator with A- 
hole excitations. These modifications to the vector me- 
son masses and pion propagator would traditionally be 
included in the chiral three-nucleon contact interaction 
and the 3NF due to intermediate A states, respectively. 

Using realistic NN interactions in many-body pertur- 
bation theory is problematic due to the strong short dis- 
tance repulsion in relative S states. The modern solu- 
tion is to integrate out the high momentum components 
of the interaction in such a way that the low energy 
physics is preserved. The details for constructing such 
a low momentum interaction, Viow-fc, are described in 
[i^ls^l- We define Mow-k through the T-matrix equiva- 
lence T(p',p,p2) = Tiow-k(p',P,P^) for {p',p) < A, where 
T is given by the full-space equation T = Vnn + Vn nqT 
and Tiow-k by the model-space (momenta < A) equation 
Tiow-k = Mow-k + Viow-kffJiow-k- Here Vnn represents 
the Bonn-B NN potential and A is the decimation mo- 
mentum beyond which the high-momentum components 
of Vnn are integrated out. Since pion production starts 
around i?iab — 300 MeV, the concept of a real NN po- 
tential is not valid beyond that energy. Consequently, 
WG choose A « 2.0 fm^^ thereby retaining only the in- 
formation from a given potential that is constrained by 
experiment. In fact for this A, the Viow-k derived from 
various NN potentials are all nearly identical [sot . 

We use the folded diagram formalism to reduce the 
full-space nuclear many-body problem H'iln = En'^n 
to a model space problem HcsXm = E„iXm as de- 
tailed in [aH. Here H = Hq + V, H^s = Hq + V^s, 
En = En{A = 14) - Eo{A = 16, core), and V denotes 
the bare NN interaction. The effective interaction Vcs 
is derived following closely the folded-diagram method 
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detailed in [34I. A main difference is that in the present 
work the irreducible vertex function (Q-box) is calculated 
from the low- momentum interaction T^ow-fci while in (3^ 
from the Brueckner reaction matrix (G-matrix). In the 
Q-box we include hole-hole irreducible diagrams of first- 
and second-order in Viow-fc- Previous studies (11,134, 35 1 
have found that Viow-k is suitable for perturbative cal- 
culations; in all of these references satisfactory converged 
results were obtained including terms only up to second 
order in Vlow-fc- 

Our calculation was carried out in _7j-coupling where 

in the basis |p3/^2 3/2^1/2 1/2} "-"^^ must diagonalize 




e 
2e 



(4) 



to obtain the ground state of ^"^N (and a similar 2x2 ma- 
trix for i-^C). We used e = E{p-^^) - E{p-J^) = 6.3 MeV, 
which is the experimental excitation energy of the first 
I state in '^^N. One can transform the wavefunctions to 
iS'-coupling, where the ^'^C and ^''N ground states are 

iP^ = xYSo) + y\'Po) 

= a\^Si)+bYPi) + c\^Di) (5) 

and the Gamow- Teller matrix element A/qt is given by 

n 



^ {^f\\a{k)T+{k)U,) = -V6 {xa - y6/V3) . (6) 



"'^^O using the charge exchange reactions ^^N(d, ^He)^^C 
and ^"'N('^He, i)^'*0, and our theoretical calculations are 
in good overall agreement. The most prominent ef- 
fect we find is a robust inhibition of the ground state 
to ground state transition for densities in the range of 
0.75 — l.OnQ. In contrast, the other transition strengths 
are more mildly influenced by the density dependence 
in BRS. In Fig. [3] we show the resulting half-life of ^^C 
calculated from the MM Bonn-B potential. 
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FIG. 2: The B{GT) values for transitions from the states 
of to the ^"^N ground state as a function of the nuclear 
density and the experimental values from [s^. Note that 
there are three experimental low lying 2"*" states compared to 
two theoretical 2"*" states in the configuration. 



Since x and y are expected to have the same sign [20| , the 
GT matrix element can vanish only if a and b have the 
same sign, which requires that the {^Si\Vcs Y'Di) ma- 
trix element furnished by the tensor force be large enough 
2l| . In Table U we show the ground state wavefunctions 
of ""^"^C and ^''N, as well as the GT matrix element, cal- 
culated with the MM Bonn-B interaction. 
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0.154 
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0.956 


-0.065 


1.0 


0.737 


0.675 


0.103 


0.273 


0.956 


0.074 



TABLE I: The coefficients of the LS-coupled wavefunctions 
defined in eq. ((Sjl and the associated GT matrix element as a 
function of the nuclear density n. 



In Fig. m we plot the resulting B{GT) 



I Mr 



2J^+TI^'^GT| 

values for transitions between the low-lying states of ^"^C 
and the ^"^N ground state for the in-medium Bonn-B 
NN interaction taken at several different densities. Re- 
cent experiments [36^ have determined the GT strengths 
from the ^''N ground state to excited states of ^^C and 
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FIG. 3: The half-life of ^"^C, as a function of the nuclear 
density, calculated from the MM Bonn-B potential. 

We emphasize that the nuclear density experienced by 
p-shell nucleons is actually close to that of nuclear matter; 
in Fig. 0] we compare twice the charge distribution of ^''N 
obtained from electron scattering experiments 



37, m 



with the radial part of the Op wavefunctions, indicat- 
ing clearly that the nuclear density for Op nucleons is 
~ O.SriQ. The first excited 0+ state of ^'^N together with 
the ground states of ^^O and ^^C form an isospin triplet. 
We have calculated the splitting in energy between this 
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FIG. 4: Twice the charge distribution of taken from [stI . 
and the fourth power of the p-shell wavefunctions. 



state and the ground state of "'^''N for a range of nuclear 
densities. Our results are presented in Fig. [5l where the 
experimental value is 2.31 MeV. 



> 

















0* 














1* 














0.0 


0.25 


0.5 


0.75 


1.0 


Exp 





FIG. 5: The splitting between the Ij^ and 0^ levels in ^''N for 
different values of the nuclear density. Also included is the 
experimental value. 

In summary, we have shown that by incorporating 
hadronic medium modifications into the Bonn-B poten- 
tial the decay of ^'^C is strongly suppressed at densities 
close to that experienced by valence nucleons in ^"'C. In 
a more traditional approach such medium modifications 
would be built in through 3N forces, and we suggest 
that calculations with free-space 2N interactions supple- 
mented with 3N forces should also inhibit the GT tran- 
sition. 
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